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ABSTRACT 


The contact problem for a thin elastic reinforcement bonded 
to an elastic plate is considered. The stiffening layer is 
treated as an elastic membrane and the base plate is assumed to 
be an elastic continuum. The bonding between the two materials 
is assumed to be either one of direct adhesion or through a 
thin adhesive layer which is treated as a shear spring. The 
solution for the simple case in which both the stiffener and 
the base plate are treated as membranes is also given. The 
contact stress is obtained for a series of numerical examples. 

In the direct adhesion case the contact stress becomes infinite 
at the stiffener ends with a typical square root singularity 
for the continuum model and behaving as a delta function for 
the membrane model. In the case of bonding through an adhesive 
layer the contact stress becomes finite and continuous along 
the entire contact area. 


This work was supported by the National Aeronautics and Space 
Administration under Grant NGR 39-007-011 and by The National 
Science Foundation under Grant GK-11977. 



1. INTRODUCTION 


The purpose of this paper is to investigate the distribution 
of the adhesive shear stress in elastic plates reinforced by an 
orthotropic elastic layer. The generalized plane stress analysis 
of lap joints shows that if the stress variation in thickness 
direction is not taken into account and if the Poisson's ratios 
of the two bonded plates are equal, then the load transfer from 
one plate to the other takes place along the edges of the bonded 
region only [1,2]. This means that the contact shear between 
the layers is zero everywhere except along the line boundary of 
the contact area where its magnitude is infinite. Even when the 
two Poisson's ratios are different, the magnitude of the contact 
shear in the contact area away from immediate neighborhood of 
the boundary is found to be extremely small [2]. This result is 
rather disturbing because of the fact that in reinforcing the 
plates with uni di recti onal ly strong layers or straps the load 
transfer is expected to take place through the adhesive bond and 
the extremely large adhesive shear along the edge of the bonded 
region means that the debonding would take place under relatively 
low values of the external loads. 

In practice the thickness of the reinforcing layer is 
generally very small compared to the thickness of the base plate 
and the in-plane dimensions of the layered medium. Hence in 
analyzing the problem the reinforcing layer may be treated as a 
"membrane". On the other hand for the base plate, because of 
its relatively large thickness, a similar assumption neglecting 
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the thickness variation of stresses may not be justified and the 
base plate may have to be treated as an elastic continuum. Also 
since the nature of the bonding between the layer and the base 
plate is expected to be an important factor influencing the 
behavior of the contact stresses, it would be useful to solve 
the problem under different contact conditions. 

In this paper we will consider the problem of an elastic 
plate reinforced by an orthotropic layer under the assumptions 
that, (a) the stress state in the layered medium is independent 
of the z coordinate (see the insert in Figure 1), (b) the bend- 
ing stiffness of the layer 2 is negligible, (c) the base plate 1 
is an elastic continuum, (d) the only external load acting on 
the medium is the uniaxial tension o-| x °° = o Q away from the rein- 
forcement region, and (e) the bonding between the two materials 
may be accomplished through either a direct adhesion (insert in 
Figure 1) or an adhesive layer of finite thickness h 3 (insert in 
Figure 3). The solution for the special case in which the base 
plate is also treated as a membrane will also be given and 
results will be compared with that of the continuum solution. 

The solution for the limiting case of this problem for h-j = °° 
and h^ = 0 is given in [3-8]. 


* 

This means that the dimension of the composite medium in 
z-direction is either very small or very large compared to 
the x-dimension 2a of the stiffener. In the former case 

< a z >ave. “ °- In the ,atter e lz ’ £ 2z ' e o ’ e o * - u l°o /E l 
if the plate is under uniaxial tension a-| = a Q away from the 

stiffener and e Q = 0 if the plate is pulled through fixed 

grips sufficiently close to the stiffener. 
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2 . 


THE CASE OF DIRECT ADHESION 


Consider an isotropic elastic plate of thickness h-| stiff- 
ened by an orthotropic layer of thickness h^ and width 2a (insert 

in Figure 1). Let the base plate be subjected to a uniaxial 

00 

tension a-j x = a Q away from the stiffening region. Assuming 
that the dimensions of the plate in x and z directions are suffi- 
ciently large, the strains in z-direction may be expressed as 

v, 

e lz = e 2z = E^~ * ^ ^ 

where E-j and v-| are the elastic constants of the base plate. 

Also, for this case it may be assumed that the stresses in the 
composite medium are independent of the z coordinate. Defining 
the contact stress at y = 0 by 


° lxy (x,0) = - P(x) (2) 

and treating the stiffening layer as a membrane, the equilibrium 
of the layer and the condition of symmetry give 

a a 

h 2°2x^ x ^ = f P ( t ) dt , / p(t)dt = 0 , p(x) = - P (-x) 

x - a 

(3.a-c) 

Using (1) and the stress-strain relations 
3 u 2 1 

e 2x = 3x” = ^ °2x " v 2x a 2z ^ ' 


v l a o 


1 


'2z 


M 2z 

from (3. a) we obtain 


(<, 2z - v 


3u, 1 - v a 

2 - 2x f p(t)dt + 


3x 


h 2 E 2x 


2z a 2x ^ * 
v l v 2z°o 


(4.a ,b) 


( | x | < a ) (5) 
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For the base plate solving the related field equations under 
the boundary conditions* 




f-p(x) , 

x <a 


ly(X,0) = 0 , 

a lxy (x,°) = 

(0 . 

A 

td 

A 

X 


ly( X » _h l ) = ® 5 

°lxy^ x ’" h l ) 

= 0 


(6) 


and adding a homogeneous solution due to cr , for the displacement 
derivative at y = 0 we find 


3 

3 x 


u-j (x,0) = -jr 5 - - J - k(x,t)]p(t)dt 

i 1 ~ a 


where the bounded function k(x,t) is given by 

0 ? -2hi a 

00 4 a h, + 2 - 4 ah, - 2 e 1 

k(x.t) = / 1 1 


0 4 a 2 h 1 2 + 2 


^2aiT, 2^7 sina(t-x)da 

- e 


( 7 ) 


(8) 


Now, if the bonding along the interface y = 0 is one of direct 
adhesion, using ( 5 ) and ( 7 ) from the continuity condition 
u 2 (x) = u-|(x,0), (-a<x<a) we obtain the following singular 

integral equation to determine the unknown function p(x): 

3 3 £1 

? / dt ' 1 / k(x,t)p(t)dt + \ J p ( t ) dt 

“9 “9 X 


1 ' v l v 2z 

= — 0 t (-a<x<a ) , 

2 ( 1 - v -j 2 ) 0 

where k(x,t) is given by (8) and 
, _ a p l (1 ' v 2x v 2z> 

- h 2 e 2 x u-v ■ 


( 9 ) 


( 10 ) 


At the end points x = + a, p(x) must have integrable 


*Using, for example, a technique similar to that outlined in [ 9 ]. 
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singularities. Hence, referring to [10], it may be shown that 
the index of the integral equation is k = +1 and the solution 
contains an arbitrary constant which is determined from the 
equilibrium condition ( 3 - b ) - 


Let us now assume that h-|/a too is sufficiently small so 
that the base plate can also be treated as an elastic membrane. 
In this case, in addition to (1), (3) and (5), from the equilib- 
rium of the base plate and from the stress-strain relations we 
obtain 

a lx = °o - F7 l p(t)dt 1 

p 2 

3u, 1-v, v, a 

a ]x + , ( I x | < a ) (11. a, b) 

Substituting from (5) and (11), the continuity condition 
3u, 3up 

W- = W ■ (-a<x<a), gives 


B - A / p(t)dt = 0 , 
x 


1 ' v l v 2z 

6 ■ E, °o ’ 


, _ ,_v l , 1 ‘ v 2x v 2z ,,, . 

A ur ur • ( 1 2 . a-c ) 


h ,E, 


h 2 E 2x 


From (12) and (3.b,c) it is easily seen that 


p(x) = | [6 ( x-a ) - 6 ( x+a ) ] , (13) 

that is, p ( x ) is zero everywhere except at x = + a where it is 
infinite. For this simple model, the stresses in the stiffener 
and in the base plate are found to be 
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( | x | < a ) , 


a 2 x ( x ) = 


W. 


°lx ( 


X) =( a ° 

< a_ - 


_B 

o Ah 


1 


( i x | >a ) , 
( | x j < a ) . 


(14. a, b) 


3. ADHESION THROUGH AN ELASTIC SHEAR LAYER 

If the bonding between the base plate and the stiffener is 
accomplished by using an adhesive material with elastic proper- 
ties different than that of the plate and the stiffener, the 
solution given in the previous section will not be valid. In 
this section the effect of the adhesive layer will be taken into 
account by assuming that due to its relatively very small thick- 
ness the adhesive may be treated as a "shear spring". That is, 
referring to the insert in Figure 3 if u 2 (x) and u-|(x,y) are the 
x-components of the displacement vectors in the stiffener and in 
the base plate, respectively, and if h^ is the thickness and 
is the shear modulus of the adhesive, the continuity condition 
along the interface may be expressed by 

ho 


u -j ( x , 0 ) - u 2 ( x ) = — p ( x ) 

3 


(15) 


where p(x) is the contact shear defined by (2), and the stiffener 
is again treated as a membrane. Noting that u 2 (0) = 0, from (5) 
we obtain 

a 

ds 
o 


V -i V 9 0 1 " X 

U 9 (x) = - V 2z - ° - X + h 2x 2z - / ds / p(t)dt 


1 


h 2 E 2x 


V ’ V 2z a x + ’ ~ V 2 * v 2z Px 

E n °o x h 0 E, KX 


1 


2 2x 


l ~ ^2 y V p - X 

/ ( x- t ) p ( t ) d t , 


h 2 E 2x 


(16) 
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where P is a constant defined by 


P = / p(t)dt . 
o 


( 17 ) 


Similarly, expressing (7) in the range (0,a) and integrating, 
we find 

1 - V, a 


u-i (x ,0) = f 5 - x - — ~ - / [1 og | t+x | - 1 og [ t-x 
1 L -] "i-h « 


1 o 


- k-j (x,t)]p(t)dt 


(0<x<a) , 


(18) 


where 


o n “ 2 01 h 1 

00 4a h, + 2 - 4ah, - 2e 

k^x.i) - J — J^irr t«sa(t-x) 


^ O O -tUrlli LUII-I 

0 / „ 2 l 2 , o - 1 - 1 


a(4a h-| + 2 - e 

- cosa(t+x)]da . (19) 

Substituting now from (16) and (18) into (15) we obtain the 
following integral equation to determine the unknown function 
p(x) : 

a , x y n h 


~ /[l°g||rf| - k 1 (x,t)]p(t)dt - J / (x-t) p ( t ) d t + y _ l ) p(x) 
o o 3 1 

( 0<x<a ) (20) 


= [ — - V | - Z - ^ ]X , 


9 M 2 , o a 
2(l-v 1 ) 


where the constants X and P, and the kernel k^(x,t) are given by 
(10), (17), and (19), res pecti vely . 


If h i / a is sufficiently small for this problem too a very 
simple solution may be obtained by treating the base plate as 
well as the stiffening layer as a membrane. In this case the 
equations (3), (5) and (11) remain valid. Thus, differentiating 
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the continuity condition (15) and substituting from (5) and (11) 
we obtain the following simple i ntegrodi ff erenti.al equation to 
determine the unknown function p(x): 

B - A / p(t)dt = ^ ^ p(x) , ( 0<x<a ) (21) 

x 3 


where B and A are given by (12.b,c). 
dit ion that p(0) = 0, we obtain 


p(x) 


£ shax 
a chaa 5 


a 



Solving (21) under the con 



( 22 . a-c ) 


It may be seen that (22) reduces to (13) as h 3 ->-0. To do 
this note that for small values of h 3 we have 

£-0(h 3 -V. a - 0(h 3 ^) , §£S| = 0(e- a < a -*>) . (23) 

Thus, since lim z"^e~ ez -> 0 as z + °° for any e>0, from (22) we 
obtai n 

f 0 , (0<x<a) , 

lim p ( x ) = < 

h 3 ->-0 oo , (x = a) , 

a R 

. 1 im / p(x)dx = j , (24. a ,b) 

h 3 +0 o rt 

which is the result found in the previous section. 


4. SOLUTION OF THE INTEGRAL EQUATIONS AND NUMERICAL RESULTS 

Referring to [10], the solution of the singular integral 
equation (9) is of the form 

p ( x ) = f(x)(a 2 -x 2 ) 2 , (25) 

where f(x) is a bounded odd function. The integral equation is 
solved by using the technique described in [11]. The major 
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numerical problem here is the evaluation of the kernel k(x,t) 


given by (8) for which a modification of Filon's integration 
formula given in [12] is used. In this problem the physical 
parameters which may be varied are the thickness ratio h-j/a for 
the base plate and the dimensionless constant X defined by (10). 
In the numerical examples given in this section the following 
material constants are used: 

Base plate: Aluminum (E^ = 10 7 psi, v-j = 0.3), 

Stiffener: Boron-Epoxy Composite (E 2x = 32.4xlo 7 psi, 

v 2x =0.23, E 2z = 3.5*10^ psi, v 2z = 0.025, 
y 2 = 1 .23xl0 6 psi) , 

Adhesive: Epoxy (u 3 = 1.65xl0 5 psi, v 3 = 0.35). 


For a fixed value of h^/a = 0.25 the contact stress p(x) 
obtained from (9) for various values of X is shown in Table 1. 
Similar results for a fixed value of X = 0.5 and various values 
of h i / a are shown in Figure 1. At the end points x = + a the 
contact stress has a conventional square root singularity the 
strength of which may be character!’ zed by a constant defined by 
(see (25)) 

K = 1 i m /2{ a-x) p(x)=f(a)//a. (26) 

x-*a 

Figure 2 shows the variation of the constant K with h^/a and X. 
The values of K corresponding to the results given in Table 1 
are shown in Table 2 where 


I V -I V 2 7 

K = '- 4 1 

o , 2x o 
2(l-v 1 ) 
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Tabl e 1 . 


Distribution of contact stress for a stiffener 
directly bonded to an elastic plate (Materials: 

1 Aluminum, 2 Boron-Epoxy composite, h-j/a = 0.2 


2(l-v 1 Z ) p ( t i ) 



9 . 9862953E-01 
9 . 8768834E-01 
9 . 6592583E-01 
9 . 3358043E-01 
8.91 00652E-01 
8. 3867057E-01 
7.771 4596E-01 
7.071 0678E-01 
6.2932039E-01 
5 . 4463904E-01 
4 . 5399050E-01 
3 . 5836795E-01 
2 . 5881 905E- 01 
1 . 5643447E-01 
5 . 2335956E-02 


4 . 1 003254E+00 
1 . 1 6027 31 E + 00 
5 . 3547896E-01 
2. 6658982E-01 
1 . 2799551 E- 01 
5.4958676E-02 
1 . 95281 45E-02 
5.1 368259E-03 
7.30091 77E-04 
9.8345780E-05 
9 .41 63598E-05 
2 . 9878907E-05 
5.4920568E-06 
4 . 0977364E-07 
6 . 06851 20E-08 


4 . 7358320E+00 
1.4171 61 5E+00 
7 . 00671 46E-01 
3 . 7 500242E-01 
1 .94031 52E-01 
9 . 05031 1 6E-02 
3 . 5737633E-02 
1 . 1 1 23089E-02 
2 . 4095920E- 03 
2 . 0800832E-04 
8.611 4432E-05 
4 . 7722433E-05 
1 . 3046450E-05 
2.25741 34 E- 06 
2.0693287E-07 


5 . 01 47577E+00 
1 . 531 2509E + 00 
7 . 76601 1 2E-01 
4.271 5300E-01 
2 . 274801 8E-01 
1 . 09581 86E-01 
4. 5058720E-02 
1 . 4889820E- 02 
3.61 96458E+03 
4 . 9638807E-04 
4 . 5557425E-05 
5.02501 21 E- 05 
1 .6648501 E- 05 
3.51 99361 E- 06 
4 . 5071 1 92E-07 


t.j/a 


9 . 9862953E-01 
9 . 87 68834 E- 01 
9 . 6592583E-01 
9. 3358043E-01 
8.91 00652E-01 
8. 3867057E-01 
7.771 4596E-01 
7.071 0678E-01 
6 . 2932039E-01 
5 . 4463904E-01 
4 . 5399050E-01 
3 . 5836795E-01 
2.5881 905E-01 
1 . 5643447E-01 
5 . 2335956E-02 



5. 1 721 770E+00 
1 . 5959247E+00 
8.20241 54E-01 
4.5769981 E- 01 
2.4751 027E-01 
1 . 21 29973E-01 
5 . 0959695E-02 
1 . 7369433E-02 
4.461 81 09E-03 
7.1678844E-04 
6.351 2786E-06 
4 . 8336798E-05 
1 .83387 00E-05 
4 . 3008397 E- 06 
6 . 2258680E-07 


5.2277471 E+00 
1 .61 87991 E + 00 
8 . 3577334E-01 
4 . 6866637E-01 
2 . 5477 51 6E-01 
1 . 2560020E-01 
5 . 31 561 64E-02 
1 . 8309233E-02 
4.7891450E-03 
8. 0595704E-04 
1 . 091 3307E-05 
■4 . 6879506E-05 
■1 .8822525E-05 
•4 . 578201 3E-06 
■6.881 1 183E-07 


5 . 2734079E+00 
1 .63761 1 0E+00 
8.485831 0 E - 0 1 
4.7774724E-01 
2.6081 9 4 3 E - 0 1 
1 .291 9782E-01 
5 . 5005804E-02 
1 . 91 07330E-02 
5.070381 0E-03 
8 .83971 27E-04 
2.6573531 E-05 
■ 4 . 5328367E-05 
■1 .91 62521 E-05 
■4.80391 58E-06 
■7.4347031 E-07 













Table 2. 


Strength of stress singularity K 
corresponding to Table 1 (— ^ = 0.25) 

a 


n 

5 

2 

1 

0.5 

1/3 

0.2 

■M 

0.2187 ' 

0.2511 

0.2653 

0.2733 

0.2762 

0.2785 


For the stiffened plate problem in which a separate adhesive 
layer is used to join the stiffener to the base plate the numer- 
ical results obtained from the solution of (20) are shown in 
Figures 3 and 4 . Here the kernel of the integral equation is 
square integrable. Hence the equation is that of a Fredholm 
integral equation of the second kind. Consequently, its solution 
p ( x ) will be bounded and continuous in the closed interval 0<x£a . 
Because of the logarithmic term in the kernel, from the integral 
equation (20) .it can be shown that at the end point x = a 
dp/dx will be infinite. Figure 3 shows the effect of h-j/a on 
p(x) for a constant X = 5 (which, for the materials under con- 
sideration, corresponds to h^/a = 0.0336) and h^/a = 0.004. For 
h-j/a = 0.25 the effect of X is seen in Figure 4. The figures 
indicate that there is a severe stress concentration at the end 
point x = a and for relatively small h-j/a, because of the "bend- 
ing" of the base plate, there is a sign reversal in the shear 
stress p ( x ) acting through the adhesive. The same stress rever- 
sal effect (in a smaller degree) is also observed in the direct 
adhesion (i.e., h^ = 0) case shown in Table 1. 

Figure 5 shows the comparison of the results obtained by 
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using three different models. Curve (a) corresponds to the 
direct adhesion (h 3 = 0) problem in which the stiffener is 
treated as a membrane and the base plate is an elastic continuum 
(equation (9)). Curve (b) contains the effect of the adhesive 
layer in addition to the assumptions made for (a) (equation (20)). 
Curve (c) is obtained from (22) where both the stiffener and the 
base plate are treated as membranes and the bonding is through 
an adhesive layer. The fourth model (i.e., two membranes bonded 
without any adhesive) would simply give a peak of infinite magni- 
tude at x = a and zero values elsewhere (equation (13)). The 
numerical solutions of the integral equations (9) and (20) indi- 
cate that as h -j / a and f^/a decrease the convergence becomes 
slower and hence computations become costlier. On the other 
hand Figure 5 shows that even for a relatively large h-j/a the 
difference between the membrane solution (c) and continuum base 
plate solution (b) may not be considered high enough to justify 
the elaborate and complicated analysis required by the latter. 

For smaller values of h^/a ( i = 1 , 2 ) the difference would, of 
course, be even smaller. Hence, it may be concluded that in the 
type of problems discussed in this paper the closed form solution 

given by (22) may give adequate results provided h-j/a and h 2 /a 

* 

are "sufficiently small" . The results also show that in the 


*The results given in this paper may be a useful guide in attempt- 
ing to interpret the phrase "sufficiently small" quantitatively. 
For example, from these results it is clear that for the practi- 
cal problems considered in [13], the membrane assumption would 
be perfectly adequate. Another model for problems of this kind 
would be the treatment of the elastic layers as "plates" with 
certain bending stiffnesses and the approximation of the adhesive 
by a combination of uncoupled shear and tension springs. This 
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presence of an adhesive layer the direct adhesion model de- 
scribed in section 2 is totally unsuitable for the analysis 
of bond rupture. 


would give rise to normal as well as shear contact stresses 
along the contact area. For geometries such as lap joints the 
concentration of the normal stresses would be quite signifi- 
cant, playing an important role in bond rupture studies. For 
the symmetric geometries, however, one would expect that the 
dominant contact stress would be the shear stress which would 
be approximately the same as that obtained from the membrane 
theory . 
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Figure 1. Distribution of contact stress for the direct 
adhesion case. Materials: 1 Aluminum, 2 Boron-Epoxy Composite* 

1 - v,v 0 _ 


0 


2.0 1/A 


4.0 


Figure 2. Strength of stress singularity for the direct 

1 - v,v 2 

adhesion case, = — o fi. 

0 2(l-v, 2 ) 0 




Figure 3. Distribution of contact stress for the case of 
bonding through an adhesive layer. Materials: 1 Aluminum 
2 Boron-Epoxy Composite, 3 Epoxy; X = 5, h 3 /a = 0.004, 

P 0 - o 0 (l -v,v 2z )/(2(1-v, 2 )). 



Figure 4. Same as Figure 3 for constant h-j/a = 0.25 and 
h Q /a = 0.004 and varying A; p Q = a Q ( 1 - v 2z )/ ( 2( 1 - V] 2 ) ) 




Figure 5. Contact stress distribution calculated from various 
models: (a) Stiffener: membrane. Base Plate: elastic continuum. 

Contact: direct adhesion; (b) Stiffener: membrane, Base Plate: 
elastic continuum. Contact: bonding through an adhesive layer; 
(c) Stiffener and Base Plate: membrane. Contact: bonding through 
an adhesive layer. A = 5, a/h^ = 4, h /a = 0.004. 


/ 



